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Abstract. The phase space representation for a g-deformed model of the quantum 
harmonic oscillator is constructed. We have found explicit expressions for both the 
Wigner and Husimi distribution functions for the stationary states of the q-oscillator 
model under consideration. The Wigner function is expressed as a basic hypergeometric 
series, related to the Al-Salam-Chihara polynomials. It is shown that, in the limit 
case h ^ {q ^ 1), both the Wigner and Husimi distribution functions reduce 
correctly to their well-known non-relativistic analogues. Surprisingly, examination of 
both distribution functions in the g-deformed model shows that, when q <^ 1, their 
behaviour in the phase space is similar to the ground state of the ordinary quantum 
oscillator, but with a displacement towards negative values of the momentum. We 
have also computed the mean values of the position and momentum using the Wigner 
function. Unlike the ordinary case, the mean value of the momentum is not zero and 
it depends on q and n. The ground-state like behaviour of the distribution functions 
for excited states in the q-deformed model opens quite new perspectives for further 
experimental measurements of quantum systems in the phase space. 

PACS numbers: 03.65.-w; 03.65.Ta; 42.50.Dv; 02.30.Gp 

Submitted to: J. Phys. A: Math. Gen. 
1. Introduction 

The harmonic oscillator concept occupies a central position in science and engineering 
due to its simplicity and exact solubihty in both classical and quantum descriptions. 
Today it appears in mechanics, electromagnetism, electronics, optics, acoustics, 
astronomy, nuclear theory, to name just a few [1]. 

In the quantum approach the harmonic oscillator is one of the exactly solvable 
problems studied in detail due to its considerable physical interest and applicability [2]. 
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To examine the applicability of a system modelled by a harmonic oscillator, one wishes 
to compare theoretical results about position and momentum with measurements. In 
practical quantum applications however, one is faced with measurability problems of 
the quantum system states due to the uncertainty principle of quantum physics. For 
this reason, it is common to work with area elements in phase space, whose size is 
not smaller than Planck's constant, and obtain an expression for the joint probability 
distribution function of momentum and position. Then, having a certain method 
for the measurement of this quantity, the gathered results will provide a complete 
characterization of the quantum system under consideration, and allow the use of more 
transparent classical language to examine its properties. 

The Wigner distribution function is the main theoretical tool to construct the phase 
space of a quantum system. It is the closest quantum mechanical analogue of the 
classical probability distribution over the phase space, i.e. limW {p,x) = p{p,x). This 
distribution function was first proposed by Wigner to study the quantum correction for 
a thermodynamic equilibrium [3]. Quantum harmonic oscillators are one of first systems 
where analytical expressions for the Wigner distribution function were calculated, 
for stationary states as well as for states of the thermodynamic equilibrium. These 
expressions are important due to experimental measurements done in the past two 
decades [4,5]. 

For example, the novel experimental method of optical homodyne tomography 
allows the measurement of the Wigner function for squeezed, vacuum, one- and two- 
photon states [4,6,7]. Note also experiments where measurements of the Wigner function 
of Fock states were performed on vibrational states of a trapped Be^ ion [8]. In 
other words, advanced measurement techniques for the Wigner function changed its 
status from an "only theoretically computable expression" to that of an "also directly 
measurable quantity". 

All of this, however, is still not applicable to generalizations of the quantum 
harmonic oscillator, in particular the so-called g-deformed harmonic oscillator. The 
origin of g-deformed oscillators can be traced back to the work of Iwata [9], who 
generalized the Heisenberg commutation relation and solved the relevant eigenvalue 
problem. Later on, in the 70's g-deformed generalizations of the harmonic oscillator were 
quite popular and since the end of the 80's they are back into fashion. It is interesting to 
observe that g-deformed harmonic oscillators models are also simple and exactly solvable. 
Moreover, the appearance of an extra parameter q gives rise to additional application 
opportunities compared to the non-relativistic harmonic oscillator. Deformed oscillators 
have found applications in the theory of quons [10], statistical physics [11], generalized 
thermodynamics [12] and in models describing a small violation of the Pauli exclusion 
principle [13], to name a few. However, in spite of these applications, many aspects of 
the g-oscillator model have not been developed. One of the main reasons restricting its 
wide applications is the absence of an exact expression in terms of phase space. In this 
context, it is necessary to note the work in [14-16], where the possibility of obtaining 
the Wigner function associated with the g- Heisenberg commutation relation is discussed. 
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However, in spite of attempts to derive an analytical expression of the Wigner function 
for the g-deformed harmonic oscillator models, up till now this goal was not achieved. 
The main problem is related to difficulties for the right definition of the momentum and 
position operators and to transformation formulae between them in the g-case. 

If one imagines modern quantum physics as a puzzle, then due to the absence of 
an explicit description of g-deformed harmonic oscillators in phase space, some pieces 
of the puzzle are missing. In this paper, we examine some missing pieces of this puzzle 
and present explicit expressions for the Wigner and Husimi distribution functions for 
the stationary states of a g-deformed harmonic oscillator model. We also analyse these 
expressions for various values of the parameter q and prove that they reduce to the 
well-known expressions for the non-relativistic harmonic oscillator in the limit q ^ 1. 

Our paper is structured as follows: in Section 2, we provide basic information 
about the definition of the Wigner function and its Gauss smoothed analogue. We 
present their explicit expressions for the stationary states of the non-relativistic linear 
harmonic oscillator. Section 3 is devoted to a model of the g-deformed oscillator, whose 
wave functions are expressed in coordinate and momentum spaces by Rogers-Szego and 
Stieltjes-Wigert polynomials respectively. We present exphcit expressions for the Wigner 
and Husimi distribution functions for the stationary states of the g-deformed oscillator 
in Section 4. Discussions and Conclusions are given in Section 5 and Section 6. 

To end this introduction, we collect some common notation and known results for q- 
series. Many of the expressions encountered in this article will be expressible in terms of 
basic hypcrgcomctric scries (g-scries), for which we use the standard notation of Gasper 
and Rahman [17]. The g-shifted factorial is defined as: 

n-l 

(a;g)o = l, and (a; g)„ = - a g*"), for n = 1, 2, . . . , oo, 

fc=0 

and (oi, . . . , a^; g)„ = (oi; g)„ • • • (a^; g)„. A basic hypergeometric series is defined as 

,(f)J «i'«2,---,ar \ _Y (ai,a2,---,a.;g)fc V^'"' 
\^ 61,62, '' y ^ (g,6i,62,...,6s;g)fc ^ ^ 

In most cases, we will have that r — s + 1. Moreover, a basic hypergeometric series 
is terminating if one of the parameters aj {j = 1, . . . ,r) equals g~" with n a nonnegative 
integer. The single most important summation formula for basic hypergeometric series 
is the g-binomial theorem: 

In the terminating case, i.e. when a — g~", this reduces to 

ifol ; ^ = —7 ^ = [q z; q)n, 

\ - J {z;q)oo 

and there are no longer any convergence conditions. 
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It will also prove convenient to have a notation for the g-binomial coefficient: 



q {<i;(l)k{q;q)n-k 
and for the g-numbers: 



a 



q^l. 



2. The Wigner distribution function 

The Wigner function for stationary states of a quantum system can be obtained from 
its wave functions ipn{x) or il^n{p) (in the position or momentum representation) by the 
well-known definition [18] 



Wn{p,x) 



Wnip,x) 



2'Kh 
1 

2t^ 



— oo 
oo 



C {p-Ip') i'n {p+^p') e"P'/^dp'. 



(2.1a) 



(2.U) 



It allows one to calculate the quantum average of a physical quantity / by the 
formula 



fn 




f{p,x)Wn{p,x)dpdx, 



(2.2) 



where f{p, x) is the Weyl symbol of the operator f{p, x) [19]. 

The following explicit expression of the Wigner function for the stationary states 
of the non-relativistic linear harmonic oscillator is well-known [19,20]: 



W^'' {p,x) 



TT 



h 



exp 



huj \ 2m 



1 1 
moo X 



■ K (- 

^ \kuj \2m 



1 1 
muo X 



(2.3) 



where L„ are the Laguerre polynomials [21, Section 1.11]. 

The Wigner distribution function determined in such a way takes both negative 
and positive values, therefore it is only a quasiprobability distribution function of p 
and X. It is bounded by the restriction |W„(p, x)| < (tt/I)""'^ [22]. For this reason 
a class of Gaussian smoothed distribution functions with nonnegative behaviour was 
introduced [23-25]: 



Wn {P, X) = — 
TXh 




exp 



P 



i2 



X 



i2 



— oo — oo 



2A2 2A2 



Wn {p + p\x + x') dp'dx' > 0, 



:2.4) 



where ApA^. is a finite region of the phase plane. The simplest form of the Gauss 
smoothed Wigner function is when ApA^. = h/2, where (12.41) will be a so-called Husimi 
distribution function [26]: 



Wn (p, X) 



{27if^ hA, 



i)n {x') ■ exp 



ipx 



{X — X 

4A2 



dx' 



(2.5) 
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Equation (12.51) has the property that Wn {p, x) is restricted by the condition 
0<Wn{p,x) < {rth)-\ 

The calculation of the distribution function of the non-relativistic harmonic 
oscillator with this formula, when the parameter is taken equal to h/2muj (simplest 
case), leads to the expression [19]: 



W" {p,x) = {27ihn\y 



lip muj 
hu \ 2m 2 



■ exp 



lip muj x" 
hjjj \ 2m 2 



.(2.6) 



3. The g-deformed harmonic oscillator model 



The g-deformed harmonic oscillator model considered in this article was already 
developed in a number of papers [27-35]. It has the following creation and annihilation 
operators: 



VI - 9 

with A = Herein, /i is a deformation parameter related to a finite-difference method 
with respect to x, and 

In the rest of this article, q will always have the value e~^^^ . 

Note that an operator of the form exp(a9a;), where a is an arbitrary complex 
number, acts as follows on functions of x: 

exp{adx)f{x) = f{x + a). 

In the limit /i — > (or g — ^ 1) the operators (13.11) reduce to the well-known non- 
relativistic harmonic oscillator creation and annihilation operators [34]: 

h^o 2VA dx 

It is easily verified that 

[r,6+], = b-b+ -qb+b- = 1. 

The Hamiltonian operator H of the model is then 

H = hu;[v^bn-+[^^). 

The wave functions for this model in the ^-representation [36] are related to the 
Rogers-Szego polynomials Hn{x \ q) [37, Chapter 17]. More in particular, one has 

7/>f o (x) = CnUn {-e-''^''^ I q) e-^^\ (3.2) 

where the normalization constant c„ is given by: 



2A^ 



TT 



{-irq-'\q-.q)-^l\ 
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One can see that in the hmit for h—^0 one recovers the wave functions for the 
non-relativistic quantum mechanical harmonic oscillator in the ^-representation: 



(x) '^-Vr {^) = I x) . e-^^\ (3.3) 



The above limit can be verified by employing the following known limit relation 
between the Rogers-Szego and Hermite polynomials [34,36]: 



hm l^-y ^ j Un (-e-^-J'lg) = {y) , q = e'^"', a > 0. (3.4) 

The wave functions in the momentum and space representation are related through 
a Fourier transform: 

^f^(p) = ^ rC^«(x)e-^da;. 
y/ Ann J-oo 

Those in the p-representation are expressed through the Stieltjes-Wigert 
polynomials S'„ (x; g): 



^r"(p) = -""^Sn [q^'^'e-^; q) e-f^ , (3.5) 

where we have used the notation and normalization of [21, Section 3.27] for the Stieltjes- 
Wigert polynomials S'„ (x; q). 

Also here, one recovers the correct expression of the wave functions for the non- 
relativistic quantum mechanical harmonic oscillator in the p-representation when h ^ 0: 

(p) '-^"V^r (P) = , (p/vWn) ■ e- = , (3.6) 

which follows from the limit relation between the Stieltjes-Wigert and Hermite 
polynomials [34,36]: 

\ira(^X'\q-q)^S^{r'''e-'^y-q)=RM, q = e-^^\ a > 0. (3.7) 

Thanks to the orthogonality relations for the Rogers-Szego and Stieltjes-Wigert 
polynomials on the real axis [36], 



oo 



^ ! -Hn (-e-2-^|g) ■ {-e^"^y\q) ■ e-y'dy = (3.^ 

V'^ J q 



— oo 
oo 



— oo 



Sn ir'^'e-'"'; q) ■ Sm {r"'e~'^y- q) ■ e"^ dy = — — (3.9) 

[q\q)n<i 



both wave functions (13.21) and (13. 5p are also orthonormal: 

oo 

rn'''''{x)rJ!''{^) = / < {pWr^^ {p) dp = 5^^- (3.10) 
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Finally, let us mention that the spectrum of the Hamiltonian H is given by: 



r 

n-\ — 
2 



with n = 0, 1, 2, . . . 



4. Computation of distribution functions 

4-1. The Wigner distribution function for stationary states 

The Wigner function for stationary states of the g-oscillator can be computed 
using (|2JU|) and or fl2n and ([33]). Substitution of in (I2JU|) leads to the 
following integral: 



2nh 



JXh{2x-x') 



Here, the Rogers-Szego polynomials have the following explicit expression: 



" r H / \ 
nn{x\q) = 22 f. (^) ' 



< g < 1. 

Using this expression, one writes the Wigner function (14. ip explicitly as 

1 " 



(4.2) 



\k+s^-(k+s)/2 



coo 

X I e 

-oo 



k,s=0 



-I 9 



-A{x-;|-)2-A(x+|-)2g2iA/i{x-;|-)fcg-2iA/i(x+|-)Sg-ipx'/;i^^/ 

Next, we use the Gaussian integral: 



^_^a2x^+aix+ao) ^ / iLg-^o+^^/^^a , when ttg > and ai, ao G C. (4.3) 




After some trivial simplifications, one then finds the following double sum expression 
for the Wigner distribution function: 



7ih (g; q)n 



k+s 



k,s=0 



n 




n 


_k_ 




s 



qy 



-ks ^—ka* —sa 



with 



hp 

a = h 2iXhx. 

n 



, (4.4) 



(4.5) 



This expression can be further reduced to a single sum by applying the g-binomial 
theorem once. Indeed, we have: 



1 



TTh{q; q)n 
1 

7r^(g; g)„ 



s=0 



ficj \ 2m ' 2 



s=0 



qy-iie (g^'e g)r 



1 „n 2_ / I mu^x^ 

^ q ftu \ 2m 2 



TT^ (g; g), 



(e-'^^g)„2^i I ^ ;lf ;g,g"e-'^ 
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Applying now the following transformation formula [17, III. 7]: 



g ^ \ _ (c/6;g)„ / q "■,b,bzq "/c 



leads to the following expression 



Wn,q{p,x) = — ^^je V ^3^2\ ^ ;g,g • 4.6 

Using the classical definition of the Al-Salam-Chihara polynomials [21, Section 3.8] 
<5n(2/;a,p|g) = — — — 3<^2 I Q ;g,g I , y = cos^, 

the Wigner function ( 14. 6p can formally be written as 

W^n.,(p,a:) = ^^-^-e-"^^-e"«H-+ ^ Ag4cos2A/ix; g^-t^, g^^e^^ |g ) 

TT^i (g;g)„ v ; 

Quite generally, Wigner functions for stationary states satisfy the following 
relation [18]: 

oo oo 

j y W„,,(p,x)-iy„,,(p,x)dxdp=^. (4.8) 



-oo — oo 



To prove this relation in the current case, one can use the double sum expression of 
the Wigner function (14.41) . After simple calculations, using the Gaussian integral (14. 3p . 
the computation reduces to checking that 

, f 7 , E (^")' ^7-# (^-^ ^)m-f: irf ^p^(g-^'; g)n = 6^^. (4.9) 

This is now readily verified since, in the left hand side, all terms except the one 
with k = k' = m = n are zero. 

From the 3(f2 expression (14.61) for the Wigner distribution function, it is easy to 
verify that in the limit for h ^ one recovers the Wigner distribution function for the 
non-relativistic quantum harmonic oscillator. 

First, using series expansion in h one sees easily that 

j.^ (g"e-";g)fc(g"e-"*;g)fc ^ ^.^ (1 - g" +^e-")(l - g^+^e"'^*) 

> ( n- n\, J- J. > 




1 - gJ+1 



KhujJ I 2m ^ 



j=0 h^O 



2 I V- 2 



k\ 
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Secondly, using a term wise limit for the 3(y92-series, one finds that: 

I . (-ir (g-";g)fc(g"e-°;g)fc(g' 
hm Wn,q[p, x) = —e v / > lim — r 



2^2 ^ ^ 



k=0 



4-2. The Husimi distribution function 

The calculation of the Husimi distribution function fl2.5p follows more or less the same 
lines as that of the calculation of the Wigner distribution function. 
First, one computes: 

^„(x')e-'P"'/'^-("-"')'/^^' dx' (4.10) 

using the expressions (13. 2p . (14. 2 p and the Gaussian integral (14. 3p . Then, one can apply 
the g-binomial theorem. Secondly, one conjugates this expression and multiplies both 
expressions, taking into account the constant factor in (12. 5p . Finally, one arrives at: 

1 (p-a/2 \ / -a* /2 \ i_( v^,m}!l^\ 

WnAP.^) = 7h^- ^^^^^ .42,„+ 2 )^ 

2'nh {q;q)n 

with a given by (14. 5p . 

Using the same kind of limit calculation as before, one sees that the limit h ^ 
(or g ^ 1) yields 



lim Wn,g{p, X 



2nhn\ 



1 2 2 

p mijj X 

— H 

huj \ 2m 2 



.,2^2 



huj \ 2m 2 



which is indeed the correct expression for the Husimi distribution function (12.60 . 



5. Discussions 



The explicit expressions obtained for the Wigner and Husimi distribution function for 
the g-deformed harmonic oscillator allow us to explore in detail the behaviour of these 
distribution functions in the phase space as q varies. First, an analysis of the ground 
state (n = 0, vacuum state) shows that the behaviour of the non-relativistic harmonic 
oscillator and its g-generalization in the phase space is the same. In other words, there 
is no impact of the parameter q in the absence of photons, and the Wigner function of 
the ground state for the g-deformed oscillator is the well-known Gaussian centered at 
the phase space point (0,0). 

For the excited quantum states the situation is rather different. In Figs. 1 and 2 
we have given a density plot for the single [n = 1) and double {n = 2) photon states 
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h=0 (q=1) h=0.6 (q=0.84) 




Figure 1. A density plot of the Wigner function of the single photon state for 
the (/-deformed harmonic oscillator, for values of /i = 0; 0.6; 1; 1.6; 2.3; 15 
{q = 1; 0.84; 0.61; 0.28; 0.07; lO'^'^) and m = w = ^ = 1. The function changes 
dramatically when the value of the parameter q decreases. Observe also a displacement 
of the Wigner function due to the non-zero average of the g-oscillator momentum. 

of the Q'-deformed harmonic oscillator in the phase space. For simplicity, we use the 
scale m — (J — h — 1. The first plot in each of the figures is the Wigner function for 
the non-relativistic quantum harmonic oscillator {h = or q — 1). As long as h < 1 
we find a phase space that is not so different from that of the non-relativistic harmonic 
oscillator, except for the occurrence of a new peak at the point a; = and p < 0. 
This value h = 1 is important, because most of the investigations of finite-difference 
generalizations of the Schrodinger equation propose to take a finite-difference step h to 
be equal to A = h/mc, i.e. the Compton wavelength of a particle with mass m (see. 
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h=0 (q=1) h=0.6 (q=0.84) 




Figure 2. A density plot of the Wigner function of the double photon state for 
the (/-deformed harmonic oscillator, for values of /i = 0; 0.6; 1; 1.6; 2.3; 15 
{q = 1; 0.84; 0.61; 0.28; 0.07; lO'^'^) and m = w = ^ = 1. The function changes 
again dramatically when the value of the parameter q decreases. The displacement of 
the Wigner function due to the non-zero average of the g-oscillator momentum is clear, 
and one can also see the impact of the momentum average dependence on n. 

for example [28,38,39]). In this case {h — 1) one has fiu) = mc^. After this point, as 
h increases [h > 1) oi q decreases, one can see a significant impact of the parameter 
on the behaviour of the g-deformed oscillator in the phase space. First one can see a 
displacement of the probability function peak towards negative values of the momentum. 
This can be understood by computing the mean values of the g-oscillator position and 
momentum. 

The mean value of the oscillator position in stationary states requires the 
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computation of the following integral: 

oo oo 

x= / xWn,q{p,x)dxdp. (5.1) 



Using the double sum expression of the Wigner function (14.41) . one can see that the 
integrand is an odd function with respect to x, and therefore (15.11) is zero. 

The mean value of the momentum for stationary states is given by a similar 
expression: 



P 



J j pWn^q{p,x)dx<lp. (5.2) 



Again, one can use the double sum expression of the Wigner function (14. 4p to 
compute this integral, leading to: 

p = —n muoh. (5.3) 

So the average of the momentum depends on the number of photons and on the 
finite-difference step h. This is the reason for the displacement of the probability function 
towards negative values of the momentum. 

The last plots in Figs. 1 and 2 show a similar picture. In other words the limit g — > 
dramatically changes the behaviour of the Wigner function. The probability distribution 
function has apparently the form of a Gaussian, and its value does not seem to depend 
on the photon number n, only its displacement in phase space depends on n. This can 
be understood analytically by computing the limit h +oo (corresponding to g — > 0). 
In fact, for n > and for fixed finite values of x and p one finds 

lim PV„,g(p,a;) =0; (5.4) 

h—r + CO 

on the other hand, for every q (or h)\ 



oo oo 

j j Wn,q{p,x)<lxdp 



1. (5.5) 



These relations indicate that lim/j^+oo VF„^q(p, x), for n > 0, actually behaves as 
Gaussian distribution function with the peak displaced towards (— oo, 0) in the (p, x)- 
plane. 

In Figs. 3 and 4, we also present plots of the Husimi distribution function for the 
g-deformed harmonic oscillator, for n = 1 and n = 2. The behaviour of the plots in this 
case is the same as for the Wigner distribution functions. Comparing the behaviour of 
the Wigner and Husimi distributions for small g- values (0 < g <C 1) indicates that both 
of them are good approximations of a displaced vacuum state, but the Wigner function 
is defined more sharply than the Husimi function, allowing a more correct determination 
of the phase space probabilities of the quantum system under consideration. In fact, for 
the Husimi function one can verify that for large /i-values, (14. lip is approximated by 
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Figure 3. A density plot of the Husimi function of tlie single photon state for 
the g-deformed harmonic oscillator, for values of ft, = 0; 0.6; 1; 1.6; 2.3; 15 
[q = 1; 0.84; 0.61; 0.28; 0.07; lO^-^^) and m w = Ti = 1. The behaviour is 
similar to that of the Wigner function. 



In a similar way, one can compute an approximation of the Wigner distribution function 
for large /i- values from (14.61) : 

1 {p+nmhui)^ mujx^ 

TTAl 

These expressions confirm the observed behaviour of the distribution functions as a 
displaced Gaussian when h is large (or when q is close to 0). 
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Figure 4. A density plot of the Husimi function of the double photon state for 
the (/-deformed harmonic oscillator, for values of ft- = 0; 0.6; 1; 1.6; 2.3; 15 
{q = 1; 0.84; 0.61; 0.28; 0.07; lO'^^) and m = u = h = 1. 



6. Conclusions 

In recent years one has seen an increased interest in the study of quantum systems 
in phase space, both analytically and experimentally. For example, nowadays the 
development of measurement technologies allows the construction of experimental setups 
to recover the Wigner distribution function of certain Fock states. The use of analytical 
expressions for the Wigner function allows for a more detailed study of quantum 
entanglement. 

We have constructed the Wigner and Husimi distribution functions for the q- 
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deformed linear harmonic oscillator, whose wave functions of the stationary states are 
expressed by means of Rogcrs-Szego and Stieltjes-Wigert polynomials in the position 
and momentum representations. The analytical expression of the Wigner distribution 
function for the stationary states can formally be expressed by means of an Al-Salam- 
Chihara polynomial. The Husimi distribution function for the stationary states is 
simpler and is expressed by means of g'-shifted factorials. Using proper limits for basic 
hypergeometric series we have shown that the Wigner and Husimi functions of the 
stationary states reduce to those for the non-relativistic quantum harmonic oscillator, 
when — > 1. 

Examination of both distribution functions in the g- model shows that, when q 
tends to 0, their behaviour in phase space is similar to the ground state of an ordinary 
quantum oscillator, but with a displacement of momentum towards negative values. We 
have computed the mean values of the position and momentum for the Wigner function 
and, unlike the ordinary case, the mean values of the momentum are not zero but depend 
on q and n. The ground-state like behaviour of the distribution functions for excited 
states in the q-case opens up new perspectives for further experimental measurements 
of quantum systems in the phase space. 
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